Relativistic second-order perturbations of the Einstein-de Sitter Universe 



Sabino Matarrese,^ Silvia MoUerach^ and Marco Bruni'^ 
^ Dipartimento di Fisica "G. Galilei", Universita di Padova, via Marzolo 8, 35131 Padova, Italy 
^ Departamento de Astronomia y Astrofisica, Universidad de Valencia, 4.6100 Burjassot, Valencia, Spain 
^SISSA - International School for Advanced Studies, Via Beirut 2-4-, 34014 Trieste, Italy 

(July 25, 1997) 



We consider the evolution of relativistic perturbations in the Einstein-de Sitter cos- 
mological model, including second-order effects. The perturbations are considered in 
two different settings: the widely used synchronous gauge and the Poisson (general- 
ized longitudinal) one. Since, in general, perturbations are gauge dependent, we start 
by considering gauge transformations at second order. Next, we give the evolution of 
perturbations in the synchronous gauge, taking into account both scalar and tensor 
modes in the initial conditions. Using the second-order gauge transformation previ- 
ously defined, we are then able to transform these perturbations to the Poisson gauge. 
The most important feature of second-order perturbation theory is mode-mixing, which 
here also means, for instance, that primordial density perturbations act as a source 
for gravitational waves, while primordial gravitational waves give rise to second-order 
density fluctuations. Possible applications of our formalism range from the study of the 
evolution of perturbations in the mildly non-linear regime to the analysis of secondary 
anisotropics of the Cosmic Microwave Background. 
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I. INTRODUCTION 



The study of the evolution of cosmological perturbations is of primary importance for understanding 
the present properties of the large-scale structure of the Universe and its origin. This study is usually 
performed with different techniques and approximations, depending on the specific range of scales under 
analysis. So, for scales well within the Hubble radius, the analysis of the gravitational instability of 
collisionless matter is usually restricted to the Newtonian approximation. As seen in the Eulerian 
picture, this approximation basically consists in adding a first-order lapse perturbation 2(^g/c^ to the 
line element of a matter-dominated Friedmann-Robertson-Walker (FRW) model, while keeping non- 
linear density and velocity perturbations around the background solution. The peculiar gravitational 
potential ipg is determined by the dimensionless matter-density contrast 5 via the cosmological Poisson 
equation, V^iy9g = AnGa'^gbS, with gi, the background matter density and a the FRW scale-factor. The 
fluid dynamics is then studied by accounting for mass and momentum conservation, to close the system 
(see, e.g., Ref. This procedure is thought to produce accurate results on scales much larger than 
the Schwarzschild radius of collapsing bodies but much smaller than the Hubble horizon, where ^gj <? 
keeps much less than unity, while the peculiar matter flows never become relativistic. The first-order 
matter perturbations obtained with this Newtonian treatment can be shown to coincide with the results 
of linear general relativistic perturbation theory in the so-called longitudinal gauge ||^. To second order, 
however, the comparison is made non-trivial by the occurrence of non-linear post-Newtonian (and higher 
order in 1/c) terms in the relativistic theory (see also Refs. Some, but not all, of the aspects of the 

relativistic treatment can be accounted for by adding an extra post-Newtonian perturbation —2ipg/c^ 
to the conformal spatial metric, an extension that leads to the so-called weak-field approximation (see, 
e.g. Ref. 1^). This improvement allows, for instance, a rather accurate treatment of photon trajectories 
in the geometry produced by matter inhomogeneities, as required in the study of gravitational Icnsing 
by cosmic structures (see, e.g., Ref. and other applications. It is worth mentioning that the full 
post-Newtonian line-element in Eulerian coordinates would also include non- vanishing shift components 
(see, e.g., Ref. 0). A second-order perturbative expansion starting from this metric would lead to the 
same result of our Poisson-gauge approach discussed below, with the obvious exception of those terms 
which are post-post- Newtonian or higher in a 1/c expansion. 

^From the point of view of the Lagrangian frame of the matter, corresponding to our synchronous and 
comoving gauge below, the Newtonian approach is quite different: the 'Newtonian Lagrangian metric' 
can be cast in a simple form, where the spatial metric tensor is written in terms of the Jacobian matrix 
connecting Lagrangian to Eulerian coordinates. According to this approach, post-Newtonian terms 
of any order appear as spatial metric perturbations over this 'Newtonian background' Q]. Without 
discussing the long list of cosmological approximation schemes which have been proposed to follow 
the non-linear dynamics of collisionless matter in the Newtonian framework, let us only mention the 
celebrated Zel'dovich approximation which is strictly related to the Lagrangian Newtonian approach. 
Various extensions of Zel'dovich theory to the relativistic case have been proposed in the literature; all 
of them however require either global or local symmetries, thereby preventing their application to the 
cosmological structure formation problem. A relativistic formulation of the Zel'dovich approximation, 
assuming no limitations on the initial conditions, is instead introduced in Ref. 

So far the list basically covers all those methods which have been devised to follow non-linear structure 
formation by gravitational instability in the Universe, with the only possible exception of a few relevant 
exact solutions of Einstein's field equations, such as the Tolman-Bondi one and some of the Bianchi 
and Szekeres solutions (see, e.g., Ref. and references therein for a review). These exact solutions, 
however, have only limited application to realistic cosmological problems. 

The study of small perturbations giving rise to large-scale temperature anisotropics of the Cosmic 
Microwave Background is instead usually treated with the full technology of first-order relativistic 
perturbation theory, either in a gauge-invariant fashion or by specifying a suitable gauge. On small and 
intermediate angular scales, however, where the description in terms of first-order perturbation theory 
is no longer accurate, second-order metric perturbations can play a non-trivial role and determine new 
contributions, such as those leading to the non- linear Rees-Sciama effect jl^. In such a case a fully 
general relativistic treatment is needed, such as that recently put forward by Pyne and Carroll and 
implemented in second-order perturbation theory by MoUerach and Matarrese ||l3| . 

The aim of the present paper is to provide a complete account of second-order cosmological perturba- 
tions in two gauges: the synchronous and comoving gauge and the so-called Poisson one, a generalization 
of the longitudinal gauge discussed by Bertschinger and Ma and Bertschinger . The former was 
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chosen here because of the advantages it presents in performing perturbative calculations. The lat- 
ter, on the other hand, being the closest to the Eulerian Newtonian picture, allows a simpler physical 
understanding of the various perturbation modes. The link between these gauges is provided by a 
second-order gauge transformation of all the geometrical and physical variables of the problem. The 
general problem of non-linear gauge transformations in a given background spacetime has been recently 
dealt with by Bruni et al. (see also and will be shortly reviewed in the following section. 

The range of applicability of our general relativistic second-order perturbative technique is that of 
small fluctuations around a FRW background, but with no extra limitations. It basically allows to 
describe perturbations down to scales which experience slight departures from a linear behavior, which, 
in present-day units, would include all scales above about 10 Mpc in any realistic scenario of structure 
formation. Accounting for second-order effects generally helps to follow the gravitational instability on 
a longer time-scale and to include new non-linear and non-local phenomena. The advantage of such a 
treatment is precisely that it enables one to treat a large variety of phenomena and scales within the 
same computational technique. 

The literature on relativistic second-order perturbation theory in a cosmological framework is not so 
vast. The pioneering work in this field is by Tomita 0], who, back in 1967, performed a synchronous- 
gauge calculation of the second-order terms produced by the mildly non-linear evolution of scalar per- 
turbations in the Einstein-de Sitter universe. Matarrese, Pantano and Saez Q obtained an equivalent 
result, but with a different technique, in comoving and synchronous coordinates. Using a tetrad formal- 



ism, Russ et al. 18 1 recently extended these calculations to include the second-order terms generated 



by the mixing of growing and decaying linear scalar modes. Salopek, Stewart and Croudace [|19|] applied 
a gradient expansion technique to the calculation of second-order metric perturbations. The inclusion 
of vector and tensor modes at the linear level, acting as further seeds for the origin of second-order 
jerturbations of any kind (scalar, vector and tensor), has been, once again, first considered by Tomita 



In this paper we study the second-order perturbations of an irrotational coUisionlcss fluid in the 
Einstein-de Sitter background, including both growing-mode scalar perturbations and gravitational 
waves at the linear level. The plan of the paper is as follows. In the next section we summarize the 
results of Ref. regarding non-linear gauge transformations for perturbations of any given background 
spacetime. In Section III we consider perturbations of a generic flat Robertson- Walker model, and we 
give the transformations between any two given gauges, up to second order. Section IV is devoted to 
the study of the evolution of perturbations in the synchronous gauge, in the specific case of irrotational 
dust in the Einstein de-Sitter background. In Section V we apply the formulas obtained in Section III 
to obtain the transformations between the synchronous [|l| and the Poisson (generalized longitudinal 
1^) gauge. Using these transformations, in Section VI the results of Section IV are transformed to 
the Poisson gauge. Section VII contains a final discussion. Appendix A reviews some mathematical 
results obtained in Ref. [ p^ and used in Section II; Appendices B and C contain useful formulas in 
the synchronous gauge, used in Section IV. Appendix D contains formulas used to obtain some of the 
Poisson-gauge results in Section VI. 



II. GAUGE DEPENDENCE AT SECOND AND HIGHER ORDER 



The idea underlying the theory of spacetime perturbations is the same that we have in any perturba- 
tive formalism: we try to find approximate solutions of some field equations, regarding them as 'small' 
deviations from a known exact background solution. The basic difference arising in general relativity, 
or in other spacetime theories, is that we have to deal with perturbations not only of fields in a given 
geometry, but of the geometry itself. 

The perturbation AT in any relevant quantity, say represented by a tensor field T, is defined as the 
difference between the value T has in the physical spacetime (the perturbed one), and the value Tq 
the same quantity has in the given (unperturbed) background spacetime. However, it is a basic fact 
of differential geometry that, in order to make the comparison of tensors meaningful at all, one has to 
consider them at the same point. Since T and Tq are defined in different spacetimes, they can thus be 
compared only after a prescription for identifying points of these spacetimes is given. A gauge choice 
is precisely this, i.e., a one-to-one correspondence (a map) between the background and the physical 
spacetime. A change of this map is then a gauge transformation, and the freedom one has in choosing it 
gives rise to an arbitrariness in the value of the perturbation of T at any given spacetime point, unless 
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T is gauge-invariant. This is the essence of the 'gauge problem', which has been discussed - mainly 
in connection with linear perturbations - in many papers p^-[2^ and review articles [ p7|j2^ , following 
different approaches. 

In order to discuss in depth higher-order perturbations and gauge transformations, and to define 
gauge invariance, one needs to formalise the above ideas, giving a precise geometrical description of 
what perturbations and gauge choices are. In a previous paper [ p5[ (see also p6[|) we have considered 
this problem in great detail, following in the main the approach used in Refs. ||2^,^,^. Instead 
of directly considering perturbations, we have first looked upon the geometry of the problem in full 
generality, taking an exact (i.e., non perturbative) point of view; after that, we have expanded all 
the geometrical quantities in appropriately defined Taylor series, thus going beyond the usual linear 
treatment. 

In this section we shall summarize the main results obtained in Ref. |l5[| . Appendix A explains in 
some more detail how to proceed in the calculations. Here and in the following Greek indices fijV, . . . 
take values from to 3, and the Latin ones . . . from 1 to 3; units are such that c—1. 

We finally remind here some basics about the Lie derivative along a vector field ^, which will be 
useful in the following. The Lie derivative of any tensor T of type {p, q) (a tensor with p contravariant 
and q covariant indices, that we omit here and in the following) is also a tensor of the same type (p, q). 
For a scalar /, a contravariant vector Z and a covariant tensor T of rank two, the expressions of the 
Lie derivative along ^ are, respectively: 

£if = ; (2.1) 

£^Zi' = Z^l^"" ; (2.2) 

Expressions for any other tensor can easily be derived from these. A second or higher Lie derivative is 
easily defined from these formulas; e.g., for a vector we have £^Z = £^{£^Z): since one clearly sees 



from (2.2) that £(^Z is itself a contravariant vector, one needs only to apply (2.2) two times to obtain 
the components of £^Z. Similarly, one derives expressions for the second Lie derivative of any tensor. 



A. Gauge transformations: an exact point of view 



A basic assumption in perturbation theory is the existence of a parametric family of solutions of 
the field equations, to which the unperturbed background spacetime belongs p9| . In cosmology and 
in many other cases in general relativity, one deals with a one- parameter family of models M\] A is 
real, and A = identifies the background A^o- On each there are tensor fields T\ representing the 
physical and geometrical quantities (e.g., the metric). The parameter A is used for Taylor expanding 
these T\; the physical spacetime M-\ can eventually be identified by A = 1. The aim of perturbation 
theory is to construct an approximated solution to tWa- 

Each one-to-one correspondence between points of A^q and points of A^a is thus a one-parameter 
function of A: we can represent two such 'point identification maps' p3[ as i/'A and ip\ (for a depiction 
of this and the following, see Fig. 1 and 2 in [^). Suppose that coordinates have been assigned on 
the background A^o, labeling the different points. A one-to-one correspondence, e.g. Va, carries these 
coordinates over A^a, and defines a choice of gauge: therefore, it is natural to call the correspondence 
itself 'a gauge'. A change in this correspondence, keeping the background coordinates fixed, is a gauge 
transformation ||2^. Thus, let p be any point in A^o, with coordinates x^{p), and let us use the gauge 
ipx'- O = ip\{p) is the point in A4a corresponding to p, to which ip\ assigns the same coordinate labels. 
However, we could as well use a different gauge, tpx, and think of O as the point of A^a corresponding to a 
different point q in the background, with coordinates 5;^: then O = 4'xip) — fxi^)- Thus, the change of 
the correspondence, i.e. the gauge transformation, may actually be seen as a one-to-one correspondence 
between different points in the background. Since we start from a point p in Mq, we carry it over to 
O = ip\{p) in A^A, and then we may come back to q in Mq with </?a^, i.e. q — (p^^{0), the overall 
gauge transformation is also a function of A, which we may denote as <I>Aj and is given by composing 
if^^ with ipx, so that we can write q — ^x{p) tp^^ltpxip))- We then have that the coordinates 
of q, x'^{X,q) = $^(a;"(p)), are one-parameter functions of those of p, a:"(p). Such a transformation, 
that in one given coordinate system moves each point to another, is often called 'an active coordinate 
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transformation', as opposed to passive ones, that change coordinate labels to each point (see Appendix 
A). 

Now, consider the tensor fields T\ on each A4\. With the gauges (p\ and ip\ we can define, in two 
different manners, a representation on A4q of each T\: we can denote these simply by T(A) and T(A), 
respectively. These are tensor fields defined on Aio in such a way that each of them has, in the related 
gauge, the same components of T\. On the other hand, r( A) an d T(A) are related by $a, which gives 
rise to a relation between their components given by Eq. ( |A13| ). Since in each gauge we now have a 
field representing T\ on Ado, at each point of the background we can compare these fields with Tq, and 
define perturbations. In the first gauge the total perturbation is Ar(A) := T(A) — Tq, and in the second 
one is AT(A) := T(A) — Tq. This non-uniqueness is the gauge dependence of the perturbations. 

It should be noted at this point that we haven't so far made any approximation: the definitions given 
above are exact, and hold in general, to any perturbative order. 



B. Gauge transformations: second order expansion 



In order to proceed and compute at the desired order of accuracy in A the effects of a gauge trans- 
formation, we need Taylor expansions. In this respect, a crucial point is that gauge choices such as 
'ijjx and ifx form one-parameters groups with respect to A, while the gauge transformations form 
a one-parameter family that in general is not a group (see Appendix A and [ p^ , p^ for more details). 
Only in linear theory the action of $a is approximated by that of the element of a one-parameter group 
of transformations. A one-parameter group of transformations is associated with a vector field ^ and 
the congruence it generates, and therefore at first order in A the effect of <I>a on the coordinates x'^{p) 



is approximated by 



-I- A^'^, whereas for a tensor T we have T 2:iT + X£^ T, as it is well known. 



However, the fact that $a does not form a group comes into play with non linearity, and one can 
show that at second order two vector fields ^(ij and ^(2) are involved, and so on. That is, the Taylor 
expansion of a one-parameter family of transformations involves, at a given order n, n vector fields 
S(fc)i k = 1 . . .n. At second order, the expansion of the transformation i''(A) = $^(a;") between the 
coordinates of any pair of points of the background mapped into one another by gives 

A2 

'5(2) 



'(A) = x'' + Ae[;) + y 



(2.4) 



This is often called an 'infinitesimal point transformation'. From this, one can always define (again, see 
Appendix A) an as sociated ordinary (passive) coordinate transformation, Eq. ( Al9| ). Substitution of 
this into Eq. ( |A13| ) gives - after properly collecting terms - the gauge transformation for a tensor T: 



f(A) = T(A) + A£^,,,r + 
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(2.5) 



If T is a tensor of type {p, q) , the components of e ach ter m i n this formula have p contravariant and q 
covariant indices, appropriately given by the rules (2.1)-(2.3). 

A simple heuristic argument that may help understanding why two vector fields are involved in 
the second order gauge transformations (and n vectors at n-th order) is the following. In practice, we 
usually consider the gauge transformation between two given gauges, e.g., in this paper, the synchronous 
and the Poisson ones. Therefore, having the conditions that fix the two gauges, the unknowns of the 
problem are the degrees of freedom that allow us to pass from one ga uge to ano ther. Then, consider 
the usual first order gauge transformation, i.e., the first order part of (2.4), (2.5): it is clear that this 
fully determines ^(i) as a field in the background. Going to second order, it should be clear by the very 
definition of a Taylor expansion that ^(j^) itself cannot depend on A, and its contribution to second order 
is built from what we already know at first order, i.e. it can only give a quadratic contribution. On 
the other hand, since at second order in both gauges we have new degrees of freedom in any quantity 
(with respect to first order), it is clear that the gauge transformation itself must contain new degrees 
of freedom: these are given by ^(2). 



Since in the two gauges we can write, respectively |30 : 

\2 

r(A) = To + XST + —S^T + 0{X^) , 

\2 

f (A) ^To + X6f + —5^f + 0{X^) , 



(2.6) 



(2.7) 
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we can substitute (2^) and (2/7) into (2^) in order to obtain, at first and second order in A, the required 
gauge transformations for 6T and S^T: 



ST = 5T + £^^^Jo, (2.8) 



S^f = S^T + 2£. ST + £'1 To + £^,,. To . (2.9) 



Eq. ( |2.8| ) is the well known result mentioned above. Eq. (2.9) gives the general gauge transformation 
for second order perturbations, and shows that this is made up of three parts: the first couples the 
first order generator of the transformation ^(i) with the first order perturbation 5T] the second part 
couples the zero order background Tq with terms quadratic in the last part couples Tq with the 
second order generator ^(2) of the transformation, in the same manner than the term £(_^j^ Tg does at 
first order, in ( |2.S| ). Eq. ( |2.9| ) also shows that there are special second order gauge transformations, 
purely due to ^(2)7 when = 0; on the other hand, a non vanishing ^(2) always affects second order 
perturbations (cf. pl| ). Finally, in some specific problems in which only effects quadratic in first order 
perturbations are important, one can consider ^(2) — 0: for instance, this is the case of back-reaction 
effects (cf. Ref. ||). 

C. Gauge invariance 

It is logically possible to establish a condition for gauge invariance to a given perturbative order n 
even without knowledge of the gauge transformation rules holdin g at t hat o rder; see p^. However, we 



shall focus here on gauge invariance to second order, using Eqs. (|2.8|), ( |2.9D . 

We need to state a clear definition of gauge invariance before giving a condition. The most natural 
definition is that a tensor T is gauge-invariant to order n if and only \i 5^T — 5^T for every k < n 
(we define 5^T := Tq, ST := S^T). Thus, a tensor T is gauge-invariant to second order if S^T = S^T 



and ST = ST. Then, from (2^) and (2_^) we see that, since ^(1) and ^(2) are arbitrary, this condition 
implies that £^To — and £^ST = 0, for every vector field ^ in the background AAq. Therefore, 
apart from trivial cases - i.e., constant scalars and combinations of Kroneker deltas with constant 
coefficients - gauge invariance to second order requires that Tg = and ST = in any gauge. This 
condition generalizes to second order the results for first order gauge-invariance that can be found in 



the literature, and is easily extended to order n; see Ref. |15 



III. GAUGE TRANSFORMATION IN A FLAT COSMOLOGY UP TO SECOND ORDER 

In view of the application that will follow in Section V, we shall now use the formulas obtained in 
the previous section to show how the perturbations on a spatially flat Robertson- Walker background 
in two different gauges are related, up to second order. This will also introduce the notation used in all 
the following sections. Here and in the following Latin indices are raised and lowered using S'^^ and Sij, 
respectively. As discussed before, we set A = 1 to describe the physical space-time. 



A. Perturbed flat Robertson— Walker universe 

We shall first consider the metric perturbations, then those in the energy density and 4-velocity of 
the matter. 

The components of a perturbed spatially flat Robertson- Walker metric can be written as 

500 = -a2(T) (^1 + 22^^^'-)^ , (3.1) 



+00 

= «^WE7^-.^^^ (3-2) 



50 

r=l 
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9ij = a^{T) 



/ + 00 \ +OC ^ 

1-2 E^^^^M ^.+E^xir , (3.3) 



\r=l 



where Xi~^^ = 0, and r is the conformal time. The functions and xtj' represent the r-th 

order perturbation of the metric. 

It is standard to use a non-local splitting of perturbations into the so-called scalar, vector and tensor 
parts, where scalar (or longitudinal) parts are those related to a scalar potential, vector parts are 
those related to transverse (divergence-free, or solenoidal) vector fields, and tensor parts to transverse 

trace-free tensors. In our case, the shift uj\ ' can be decomposed as 

,^W=5iwMII+cuW^, (3.4) 

where co^^^'^ is a solenoidal vector, i.e., d^u^^^'^ = 0. Similarly, the traceless part of the spatial metric 
can be decomposed at any order as 

where a suitable function, Xi"^^ is a solenoidal vector field, and Q'Xij^"'^ = Oj hereafter, 

Bij:=didj-^5ijV^ . (3.6) 

Now, consider the energy density g, or any other scalar that depends only on r at zero order: this 
can be written as 

-l-oo ^ 
r=l 

For the 4- velocity of matter we can write 

«'^ = ^(^o^+g^<o) • (3-8) 

In addition, u^^ is subject to the normalization condition u^^u'^g^i, = — 1; therefore at any order the time 
component t;^^^ is related to the lapse perturbation, ip^r) ■ For the first and second-order perturbations 
we obtain, in any gauge: 

^(1) = ; (3-9) 

4) = -V'(2) + 3V'(i) + 2a;f + v^'^l^^ . (3.10) 
The velocity perturbation v^^^^ can also be split into a scalar and vector (solenoidal) part: 

As we have seen in the last section, the gauge transformation is determined by the vectors ^(r). 
Splitting their time and space parts, one can write 

fir) = o^^'^ , (3.12) 

and 

^(V) = + d^"-)* , (3.13) 

with diS"-^' = 0. 
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B. First-order gauge transformations 



We begin by reviewing briefly some well-known results about first-order gauge transformations, as 
we shall need them in the following. As in Section ||, we simply denote quantities in the new gauge by 
a tilde. 



From Eq. (2.8), it follows that the first-order perturbations of the metric transform as 



(3.14) 



where g^^J is the background metric. Therefore, usin g Eq. (|2.3D , we obtain the following transformations 

(3.15) 



for the first-order quantities appearing in Eqs. (3.1)-(3.3): 

7 , / a' 

= + «(!) H — "(1) ; 



(1) =</-(!)- ^ V2/3(i) - ^ 



X. 



(1) 



where a prime denotes the d eriv at ive w ith re spec t to r. 
For a scalar from Eqs. (2.S), (3.7), and ( |2.l| ) we have 

Sg = Sg + g[Q-^a(^i) . 



For the 4- velocity w'^, we have from Eqs. (2 



Using Eqs. (2.2) and (3.5) this gives: 





-"(1) 


a 

- -"(1) -"(1 


^(1) 


= "(1) 


- P(ij - (2(1) . 



(3.16) 
(3.17) 
(3.18) 

(3.19) 
(3.20) 

(3.21) 
(3.22) 



The 4-velocity is however subject to the constraint (3.9), therefore Eq. (3.21) reduces to Eq. (3.15) 



C. Second-order gauge transformations 

We now extend these well-known transformation rules of linear metric perturbations to the second- 
order. 



Second-order perturbations of the metric transform, according to Eq. (2.9), as 



(3.23) 



This leads to the following transformations in the second-order quantities appearing in Eqs. (3.1)-( p.3| ): 
lapse perturbation 



HI) ( 2^, 



',(1) 



a 



(1)1 . a 1 



.a 
) — I 

a 





a'2\ 1 







2a 



(1) (2V'(i 



•^(1) 



(3.24) 



+^1) (a?-ep-2.W)+a^,) + ^a(,, 
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shift perturbation 



(2-1 - + ^S') + ' - 3a!^^ + ") (3.25) 



spatial metric, trace 

0(2) = 0(2) + ^(1) 
1 



2N'(i)+2-'/'(i) - - + -T "(i)--«(i) 



a 



^«(i))v%)-i(2^i)-a(\)+efi))a« 



(3.26) 



3 y^'^^y ' "^i.j ^j-i ) ^(1) 
spatial metric, traceless part 



-2 - 



~(2) _ (2) „ / (1)/ o«' (1)\ ^ , <rfe 

Xij - Xij + ^ yxij + ^~Xij I "(1) + ^Xtj^k^(i) 



+2 
+2 



(3.27) 



+2 (d[^^)^.)+D,,/3(2) 
For the energy density g, or any other scalar, we have from Eq. (2.9) 



^From this we obtain, using Eq. (2.1): 

S^g 6^g + g[o)a{2) + "(i) (£'(o)a(i) + 6'(o)"(i) + '^^s') + ^(i) {^[q)^T + "^^Q-. 

For the 4-velocity , we have from (2.9): 



(3.28) 



(3.29) 



(3.30) 



Using Eqs. ( |3.§| ) and (2.5) this gives: 



^(■2) = 4) - ^"(2) - "(2) + "(1) 



2 ^; 



^(1) - 7^(1) 



a'2 a"\ a' 



(3.31) 



«(2) = «(2) - P'(2) - <2) + 



2 w 



'(i)-7"(i)) - (^(i)-2-e('i) 



'(1) 



(2-(i). - ^(i).) - ^(1). (2-^) - e^)) + efi) (2^(1) + a'd)) ; 



(3.32) 



for the time and t he sp ace components respe ctively. Again, t he 4- v elocity is subject to u^u'^ g^v 
which gives Eq. ( |3l0| ); therefore Eq. ( pi] ) reduces to Eq. ( |3^ ). 



IV. EVOLUTION IN THE SYNCHRONOUS GAUGE 
A. General formalism 

In this section we will obtain the second-order perturbations of the Einstein-de Sitter cosmological 
model in the synchronous gauge, including scalar and tensor modes in the initial conditions. The 
synchronous gauge, that has been one of the most frequently used in cosmological perturbation theory, 
is defined by the conditions goo = — and g^t = j2l]. In this way the four degrees of freedom 
associated with the coordinate invariance of the theory are fixed. 

We start by writing the Einstein's equations for a perfect fiuid of irrotational dust in synchronous 
and comoving coordinates. The formalism outlined in this subsection is discussed in greater detail in 
Ref. . With the purpose of studying gravitational instability in the Einstein-de Sitter background, 
we first factor out the homogeneous and isotropic expansion of the universe. 

The line-element is written in the form 

ds^ = a^{T) [ - dr^ + 7y(x, T)dx'dx^] , (4.1) 

with the spatial coordinates x representing Lagrangian coordinates for the fluid elements. The scale- 
factor a(T) oc is the solution of the Friedmann equations for a perfect fluid of dust in the Einstein-de 
Sitter universe. 

By subtracting the isotropic Hubble-flow, one introduces the extrinsic curvature of constant r hy- 
persurfaces, 

= , (4.2) 

with a prime denoting differentiation with respect to the conformal time r. 

One can then write the Einstein's equations in a cosmologically convenient form. The energy con- 
straint reads 

Q 24 

^?'-^?^■^?^■ + -1? + 7^=^^, (4.3) 

where TV'j{'j) is the intrinsic curvature of constant time hypersurfaces, i.e. the conformal Ricci curvature 
of the three-space with metric 7^, and TZ = 7?.\. We also introduced the density contrast S = {g — 
8(0))/ 9(0)^ with £i(x, r) the mass density and £'(o)(t) = 3/27rGa^(r)T^ its background mean value. 
The momentum constraint reads 

= ^.3 , (4.4) 

where the vertical bar indicates a covariant derivative in the three-space with metric jij. 

Finally, after replacing the density from the energy constraint and subtracting the background con- 
tribution, the evolution equation for the extrinsic curvature reads 

+ ^d'j + + i (^t^'^^/fc - t?'^ (5^ + n\ - ^TZS'j = . (4.5) 

Also useful is the Raychaudhuri equation for the evolution of the peculiar volume expansion scalar 
1?, namely 

+ -d + ^\^\ + -^S^O . (4.6) 

An advantage of this gauge is that there are only geometric quantities in the equations, namely the 
spatial metric tensor with its time and space derivatives. The only remaining variable, the density 
contrast, can indeed be rewritten in terms of jij, by solving the continuity equation. We have 

5(x, r) = (1 + 5o(x)) [7(x, r)/7o(x)] - 1 , (4.7) 

with 7 = det 7y. We denote by a subscript without parenthesis the initial condition of the referred 
quantity. 
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B. First-order perturbations 



We are now ready to deal with the equations above at the linear level. Let us then write the conformal 
spatial metric tensor in the form 



sij ■ 



According to our general definitions we then write 



with 



Recall that at first order the tensor modes xlf^^ ^re gauge- invariant. 

As it is well known, in linear theory, scalar, vector and tensor modes are independent. The equation 
of motion for the tensor modes is obtained by linearizing the traceless part of the 1?'^- evolution equation. 
One has 



(4.8) 



(4.9) 



(4.10) 



X: 



(l)T'' 



+ -X 



(4.11) 



which is the equation for the free propagation of gravitational waves in the Einstein-de Sitter universe. 
The general solution of this equation is 



X,7^(^'^)= / '^'kexp(ik.x)xW(k,T)e^,.(k), 



(4.12) 



where e^ (k) is the polarization tensor, with a ranging over the polarization components +, x, and 
Xa\^, t) the amplitudes of the two polarization states, whose time evolution can be represented as 



X^^\k,T) = A{k)a„{k) 



3jl(fcT) 

kr 



(4.13) 



with ji the spherical Bessel function of order one and a^r (k) a zero mean random variable with auto- 
correlation function (ao-(k)ao-/(k')) = (27r) 'A:^'^(5''(k + k')^^^'- The spectrum of the gravitational wave 
background depends on the processes by which it was generated, and for example in most inflationary 
models, A{k) is nearly scale invariant and proportional to the Hubble constant during inflation. 

In the irrotational case the linear vector perturbations represent gauge modes which can be set to 
zero: xl = 

The two scalar modes are linked together via the momentum constraint, leading to the condition 



D 



"so 



CsO 



The energy constraint gives 



-Xs 



12. 



(4.14) 



(4.15) 



having consistently assumed 1- 

The evolution equation also gives an equation for the scalar modes, 

^(1)11% 1^(1)11' + W)II=_24^). 
r 3 

An equation only for the scalar mode Xs^'*" can be obtained by combining together the evolution equation 
and the energy constraint. 



(4.16) 
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■^(1)11" , 2 (i)ii'^. (1)11 



12 , . . 

= -Z^So. (4.17) 



On the other hand, by hnearizing the solution of the continuity equation, we obtain 

Si''=So-lv^{x^''-xil''), (4.18) 
which replaced in the previous expression gives 

si'y' + ^sP'--^Si'^=0. (4.19) 

This is the equation for linear density fluctuation (see, e.g., Ref. [Q), whose general solution is straight- 
forward to obtain. 

The equations above have been obtained in whole generality; one could have used instead the well- 
known residual gauge ambiguity of the synchronous coordinates (see, e.g., Refs. |0,E8[) to simplify their 

form. For instance, one could fix Xg^''" so that V^Xsq''" = — 2(5o, and thus the Xs ^svolution equation 
takes the same form as that for 6. With such a gauge fixing one obtains 

xi'^"(x,T) =x+(x)r2+x-(x)r-3 , (4.20) 

where x± set the amplitudes of the growing (-I-) and decaying (— ) modes. In what follows, we shall 
restrict ourselves to the growing mode. The effect of the decaying mode on second-order perturbations 
has been considered in Ref. [0 and in Ref. and will not be studied here. The amplitude of the 
growing mode is related to the initial peculiar gravitational potential, through x+ = where in 

turn, (f is related to Sq through the cosmological Poisson equation V^(p{x) — -%(5o(x). Therefore, 

D.,xi'^" = (^ip,.j ~ . (4.21) 

The remaining scalar mode 

0«(x,r) = ^^(x) + ^VMx) (4.22) 

immediately follows. 

The linear metric perturbation therefore reads 

7g = -yV^^.-y^..+xr- (4-23) 
With purely growing-mode initial conditions, the linear density contrast reads 

4" = ^V^^ . (4.24) 
6 



C. Second-order perturbations 



The conformal spatial metric tensor up to second order is expanded as 

7., =^.,+72 + ^7^, (4.25) 

with 

7g=-2#(5.,-t-xi?] (4.26) 

and xf = 0. 
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The technique of second-order perturbation theory is straightforward: with tlie help of the relations 
reported in Appendix B, we first substitute the expansion above in our exact fluid-dynamical equations 
(momentum and energy constraints plus evolution and Raychaudhuri equations) obtaining equations 

for 7g^j with source terms containing quadratic combinations of 7g|j plus a few more terms involving 

f2") (2) 

Sq. Next, we have to solve these equations for the modes (pi ' and Xsij terms of the initial peculiar 

gravitational potential (p and the linear tensor modes Xij^"""- 

Let us now give the equations which govern the evolution of the second-order metric perturbations. 

Raychaudhuri equation 



T T 



1 ^.,(1) 



I Tsij ^^sij 



-%k \% 



^\ isi 'sOi 



1 



energy constraint 



si "sOi 



(4.27) 



1 



/ 2 (1) , (l)'s 2 (^^'^ 1 -I- I (^^-^ 

7s I V 7sii 



, 3 (i),fe _ 1 (i)fe _ 1 

's ,k isij 2 s ,fe 'si ,j ^ 'si isj ,k 



+ - 



A I "Tsi "TsOi I 9 W ^ '^sij %0 7sOij 



(4.28) 



momentum constraint 



' 2 \ ^-J*^'* 'sifej j ' IS ,i isjk 2 ^ '-J sife 2 ^* ' v / 



evolution equation 



i(2)"+lw,(2)' \ J^ , 1 f ^(2)^" , ^J^)A + J^)'^ _ 1 (2)fe^ . , 1 (2)fc* , 1(2)^,^ _ i y^Y^")' 



_ (l)ife' (l)fc' (1)*' I 1 

7s Tsfej 2'''sfc 7sj g 



7s^ 7s/c 



(l)i/ (l)k,t. „2 (l)fc 

-7s/ 7s/ ,fc - V 7sfc 



+27 



{l)kt (l)l 



Al)r , (1) ,^ _ (1)* _ (1) A. 2^(1)'=^ (^^(1)' _ ^(1)* _ ^(1)^* 
7sj ,kt ^ 'ski ,j 1st ,jk Isij ,k j ^ ^7s ,fe I 7sj ,e 1st ,j 'sjt 



(l)fei „ (l)fci (l)t {l)kt (1) ,i (1)£ / (l)fei (l)fe,i (l)i,fe 



+^7s ,^7sjfc ^7s ,^7sj ,k + 7s jlskt + 7s£ ,fe I 7s ,. 



■7 



-7s 



s I V 7s^^ -I- 7s„ ^fe^ Z7sfc I j 7s ^ 7sm ,k Isk ,m )" j 



.7s ,m7sfe^ " 97s ,m7sfe ,^0 j + .7sfe 7s« .m" j • 



(l)fc€ ^(l)m .i 



(l)A;,m (1)^ 



(4.30) 



The next step is to solve these equations. In these calculations, we can make the simplifying assump- 
tion that the initial conditions are taken at conformal time tq = (implying also 60 =0). One can 
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start from the Raychaudhuri equation, to obtain the trace of the second-order metric tensor. (Actually, 
in order to obtain the sub-leading mode generated by linear scalar modes, we also need the energy 
constraint). The resulting expression for (j)^^ is 

= ^ (-yr'^V... + (VV)^) + ^ (^^V.. + ^^VV) + , (4.31) 

(2) (2) 

where <Pl^^y which is the part of ' generated by the presence of tensor modes at the linear level, reads 



Mi) ~ 5 



with Q(x, r) a source term whose explicit form is reported in Appendix C. 

(2) . . . (2) 

The expression for Xsij is obtained by first replacing ifig into the remaining equations and solving 
them in the following order: energy constraint — > momentum constraint — > (traceless part of the) 
evolution equation. We obtain 

xS = ^ (lQ^:i^,kj - 12¥',i,VV + 4(VV)'^i,- - y¥'''V,«^i. 

(2) f2") 

where Xs(t)y P^^* °f traceless tensor xl/j generated by the presence of tensor modes at the 

linear level and includes the effects of scalar-tensor and tensor-tensor couplings; its expression can be 
derived from the equations given in Appendix C. The transverse and traceless contribution tts^ , which 
represents the second-order tensor mode generated by scalar initial perturbations, is determined by the 
inhomogeneous wave-equation 

<y + ^<^J " V^sy = -^V^^.,, (4.34) 

with 

Sij = V^oSij + *o,y + 2 (v,i, VV - ^,ik<Pj) , (4.35) 

where 

V^^-o = -I ((VV)' - ^,^k<P'''') ■ (4.36) 

This equation can be solved using the Green method; we obtain for Wsij that 

r\ , , 4t2 
^5ii(x) + ^ 



TTijix, t) = — 5y(x) + —Tij{x) + #y (x, r), (4.37) 



where V^Tij = Sij and the remaining piece TTy, containing a term that is constant in time and another 
one that oscillates with decreasing amplitude, can be written as 

^'•i('^'^) = (^/ '^'kexp(ik.x)^5,,(k) Q - ^) , (4.38) 

with SijCk) = J (i'^xexp(— ik • x)Sij{x). 
The second-order density contrast reads 

= ^(^(VV)' + 2^'^V,i,) + ^(l5^'V,i +40^VV - 6^''^xl]^^) (4.39) 
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An important aspect of our results is that linear tensor modes (gravitational waves) can generate 
second-order perturbations of any kind (scalars, vectors and tensors). This interesting fact, which 
has been first noticed by Tomita is nicely displayed by the above formula for the mass-density 
contrast, which even in the absence of initial density fluctuations, takes a contribution from primordial 
gravitational waves. More in general, we should stress that our expressions completely determine the 
rate of growth of perturbations up to second order. 



V. FROM THE SYNCHRONOUS TO THE POISSON GAUGE 

In this section we are going to obtain the metric perturbations in the Poisson gauge by transforming 
the results obtained in the synchronous gauge in the previous section. The Poisson gauge, recently 
discussed by Bertschinger Q and Ma and Bertschinger [O, is defined by uj,(^)-'^ = = 0. Then, 

one scalar degree of freedom is eliminated from goi {u;^^'>'^= 0), and one scalar and two vector degrees 
of freedom from gij (x'-'^''" = xf^'^ = 0)- This gauge generalizes the well-known longitudinal gauge to 

(r) (r) 

include vector and tensor modes. The latter gauge, in which ljI = xlj = 0, has been widely used 
in the literature to investigate the evolution of scalar perturbations [ p8| . Since the vector and tensor 
modes are set to zero by hand, the longitudinal gauge cannot be used to study perturbations beyond 
the linear regime, because in the nonlinear case the scalar, vector, and tensor modes are dynamically 
coupled. In other words, even if one starts with purely scalar linear perturbations as initial conditions 
for the second-order theory, vector and tensor modes are dynamically generated p]. 



A. First-order transformations 



Given the perturbation of the metric in one gauge, it is easy to obtain, from Eqs. ( ^.15 )-( 3J^ ), the 



gauge transformation to the other one, hence the perturbations in the new gauge. In the particular 
case of the synchronous and Poisson gauges, we have: 

; (5.1) 

a 

a(i)=/?W'; (5.2) 

= d^" ; (5.3) 



3 



(5.4) 



(xi')" +2/3(1)) =0; (5.5) 



xrfc)+C=0; (5.6) 



Xp ij Xs ij ■ 

(5.7) 



The par amet ers a^'^\ and d^p of the gauge transformation can be obtained from Eqs . (5.E), 



(5.2), and (^.6| ) respectively, while the transformed metric perturbations follow from Eqs. ( p.l| ), (5.S 
Q, and 0. 

Once these parameters are known, the transformation rules for the energy density g or any other 
scalar, and those for the 4- velocity w^, follow trivially from Eqs. ( |3.19 ), ( 3.21[ ), and ( 3.22| ). In the 



irrotational case studied in the last section x^if = ^(i)± ^ ^ ^"^^ thus d^ = uj^^l = x''p\j = 



15 



B. Second-order transformations 



Th e mo re general transformation expressions follow straightforwardly from Eqs. (3.24)-(3.27), ( [3.29| ), 
and d^H). 



(2) 

Transforming from the synchronous to the Poisson gauge, the expression for -^p can be easily ob- 
tained from Eq. ( ^.24 ), using Eq. ( |5.2| ) and the condition d^^^ = to express all the first-order quantities 
in terms of /3*^^h 



V'^'^ = /3(i) 



+ {^T' + ^(^T^ + 2/9(1) + + 



For wi^l- and (h'^^ we get: 



,(2) 



(5.8) 



(5.9) 



(1) 



2(<^«'-H2^^^« 



-Ufl'* (oAl-^'i "'fl(l)'^ 4- 2/o(l)^,jj 2 (i)T a' 1 



a 

.a' „, 4, 



(5.10) 



For xl'^L we obtain: 



(5.11) 



+2xl;f/3^')''= + 2x«^/??''= + 2x5l)^/3^^)-'= 

-^'5,x«"/3«'"= + 2(4)^+D,;3(^)) . 

Given the metric perturbations in the synchronous gauge, these constitute a set of coupled equations 
for the second-order parameters of the transformation, a^'^\ and df'\ and the second-order metric 

perturbations in the Poisson gauge, , Wp^' , (0^ , and xl^ij • This system can be solved in the following 

(2) (2) 
way. Since in the Poisson gauge 9*Xp a — 0; we can use the fact that d^d^Xp a = and the property 



= 0, together with Eq. ( 5.11 ), to obtain an expression for V2V2/3(2)^ from which /J^^) can be 
computed: 



+4V^/3,^.^/?(;^) - ^'P^^^'P]:' + ^P^P% - §V^/3(i)V2v2/3(i) (5.12) 



'(1) 
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Then, using the condition d^x^p\j — ^^'^ substituting [3^'^^ we obtain an equation for d. 



(2). 



FinaUy, using S^Wp^^^ — and substituting /S*^^^ , we get an equation for a'^' : 



2 



,(1)' 



(5.13) 



(5.14) 



Having obtained, at least imphcitly, all the parameters of the gauge transformation t o se co nd o rder, 
one can in principle compute the metric perturbations in the Poisson gauge from Eqs. (5.8)-( 5.11 ). 

Similarly, once the parameters are known, the perturbations in any scalar and 4-vector, and in 
particular those in the energy density and in the 4- velocity of matter, follow trivially from Eqs. ( 3.29 )- 

iH). 



VI. EVOLUTION IN THE POISSON GAUGE 



We have obtained in the previous section the general gauge transformation to go from the synchronous 
to the Poisson gauge up to second order in metric perturbations. We can now apply it to the case of 
cosmological perturbations in a dust universe and compute the perturbed metric in the Poisson gauge 
from the solutions obtained in Section IV using the synchronous gauge. 



A. First-order perturbations 



For the first order, replacing Eq. ( 4.21 ) in Eq. (5^) and using Eq. (5^), we obtain that the parameters 
of the transformation are 



and d^^^* = 0, in the absence of vector modes in the initial conditions. 

For the metric perturbations we obtain from Eqs. ( |5.lD , (^.3|), ( |5.4| ) and (5.7) 

4'^ = = ^, 



(6.1) 



(1) T 

X V 



(1) 



(6.2) 

These equations show the well-known result for scalar perturbations in the longitudinal gauge and the 
gauge invariance for tensor modes at the linear level. 
The linear density contrast reads 

^2 



r(l) 



^2(^+— VV 
6 

while the first-order 4-velocity perturbation has components 

(1)0 

Vp = —(fi 



„(1)'' 



(6.3) 

(6.4) 
(6.5) 
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B. Second-order perturbations 



For the second-order parameters of the g auge tran sformation, replacing the second-order perturbed 
metric obtained in Section IV in Eqs. ( ^.12 ) to ( 5.14 ) we obtain 



.(2) - 



21 ' 



10 o 



4er 



6 V12 



with V^9o = ^'o - ^</''V,i and 



4^2 



2^ 



660 



(6.6) 



(6.7) 



where the quantities indicated by the subscript (t) stand for the contributions arising from the presence 
of tensor modes at the hnear level and are discussed in greater detail in Appendix D. 
For the perturbed metric we obtain 



(2) 



1 



a(2) _ 



10 

1 

21 



10 T 

21 



3 

3^' -860, 



(2) 
p(t) 



/,(2) 
P(t) 



(2) 

pij 



2^r 



V2, ,(2)i 



X 



(2) 

p{t)ij- 



(6.8) 



The equations determining the contribution from linear tensor modes are given in Appendix D. Note 

(2) (2) 

that the contribution to "^p and (pp from linear scalar modes can be recovered, except for the sub- 
leading time- independent terms, by taking the weak-field limit of Einstein's theory (see, e.g., Ref. ^) 
and then expanding in powers of the perturbation amplitude. 

Also interesting is the way in which the second-order tensor modes, generated by the non-linear 
growth of scalar perturbations, appear in this gauge: the transformation from the synchronous to the 
Poisson gauge has in fact dropped the Newtonian and post-Newtonian contributions, whose physical 
interpretation in terms of gravitational waves is highly non-trivial (see the discussion in Ref. [Q); what 
remains is the tensor TTy , whose evolution is governed by the equation 



40 



T 



(6.9) 



Its solution, Eq. (4.38), contains a constant term, deriving from the vanishing initia l conditions, plus a 
wave-like piece, having exactly the same form as linear tensor modes (cf. Eq. ( 4.13| )), whose amplitude 
is fixed by the source term Tlj (a quadratic combination of linear scalar modes). A more extended 
discussion of these tensor modes is given in Ref. . 

Finally, let us give the Poisson gauge expressions for the second-order density and 4-velocity pertur- 
bations. One has 



c(2) _ 



^ (^5(VV)' + 2ip'''^.j + + ^ (^-21^>,^ + 24^V 

O 

-ip'^ - 246o 



144 

T 



^*o - Q^'^xS^^ 



y(l)T»Jy(l)T 

A A.IJ 



AO Xoij 



(6.10) 



and 
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,(2)0 



„(2)» _ 



1 



10 



6' 



7 



7 



+ - o<^ 1290 - 



.,(2) 
P(t) 



6 



2t 



(6.11) 
(6.12) 



with the vectors d^^-** defined in Eq. (6.7). 

In concluding this section, let us emphasize that all the second-order Poisson gauge expressions 
obtained here are new. Only a few terms in these expressions were already known in the literature, 
based on the weak- field limit of general relativity (e.g., Ref. 0). 



VII. CONCLUSIONS 



In this paper we considered relativistic perturbations of a collisionlcss and irrotational fluid up to 
second order around the Einstein-de Sitter cosmological model. The most important phenomenon of 
second-order perturbation theory is mode mixing. An interesting consequence of this phenomenon is 
that primordial density fluctuations act as seeds for second-order gravitational waves. The specific form 
of these waves is gauge-dependent, as tensor modes are no longer gauge-invariant beyond the linear 
level. A second interesting effect is the generation of density fluctuations from primordial tensor modes. 
One can even figure out a scenario in which no scalar perturbations were initially present, but they were 
later generated, as a second-order effect, by the non-linear evolution of a primordial gravitational-wave 
background. 

The first effect, which is discussed in some detail in Ref. [Q, in the synchronous and comoving 
gauge also contains a term growing like r'' and a second one growing like r^: the first accounts for 
the Newtonian tidal induction of the environment on the non-linear evolution of fluid elements, the 
second is a post-Newtonian tensor mode induced by the growth of the shear field. The remaining 
parts of this second-order tensor mode (excluding a constant term required by the vanishing initial 
conditions) oscillate with decaying amplitude inside the horizon and describe true gravitational waves. 
Quite interesting is the fact that these are the only parts of these second-order tensor modes which 
survive to the transformation leading to the Poisson gauge. 

The second effect is less known, and was only previously considered by Tomita back in the early 70's 

One may naturally wonder whether there is any hope to detect the cosmological stochastic 
gravitational-wave background produced at second order by scalar fluctuations. It is, of course, the 
oscillating part of TTy which is relevant for earth or space detectors. The problem for these wave-like 
modes is that their energy density suffers the usual dilution caused by free-streaming inside the 
Hubble radius, while at horizon -crossing their closure density is already extremely small, ^ 
(where Sh is the rms density contrast at horizon-crossing), because of their secondary origin. More 
promising is the possibility that a non-negligible amount of gravitational radiation can be produced 
during the strongly non-linear stages of the collapse of proto-structures, an issue which would however 
require a fully non-perturbative approach. 

It should be stressed that, while many of our second-order terms had already been computed in 
the synchronous gauge, all our second-order Poisson-gauge expressions are new. This is a relevant 
result, as the latter gauge is the one which allows the easiest interpretation of the various physical 
effects. In particular, the second-order metric perturbations obtained by our method allow to compute 
self-consistently gravity-induced secondary anisotropics of the Cosmic Microwave Background. This 
calculation has been recently performed by MoUerach and Matarrese ||T^, implementing a general 
scheme introduced by Pyne and Carroll [Oj. 
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APPENDIX A: TAYLOR EXPANSION OF TENSOR FIELDS 



In this appendix we present some mathematical results used in Sec. II, concerning Taylor expansions 
of tensor fields on a manifold. These results have already been presented in [p^ , where analyticity of all 
relevant fields was assumed; they have been generalized in |^ to the case of C™ fields. The theorems 
obtained in are very general, concerning perturbation theory at an arbitrary order n. In order 

to achieve these general results it is very useful, or perhaps mandatory, to use a fully geometrical 
approach. However, for our purposes, it is useful to summarize them in terms of coordinates and tensor 
components, as we shall do in the following. We assume that all quantities are as smooth as necessary. 



1. One-parameter groups of transformations 



As disccused in Sec. II, gauge choices for perturbations entail the comparison of the tensor field 
representing a certain physical and/or geometrical quantity in the perturbed spacetime with the tensor 
field representing the same quantity in the background spacetime. Consequently, gauge transformations 
entail the comparison of tensors at different points in the background spacetime. A smallness parameter 
A is involved, so that these comparisons are always carried out at the required order of accuracy in A, 
using Taylor expansions |Q . Differential geometry tells us that the comparison of tensors is meaningful 
only when we consider them at the same point. Therefore, supposing we want to compare a tensor field 
T at points p and q, we need to define a transport law from q to p. This gives us two tensors at p, T 
itself and the transported one, which can be directly compared. 

The simplest transport law we need to consider is the Lie dragging by a vector field, which allows 
us to compare T with its pull-back T{X) (the new tensor defined by this transport). To fix ideas, let 
us first consider, on a manifold M, the comparison of tensors at first order in A (which we shall define 
shortly). Suppose a coordinate system has been given on (an open set of) M., together with a vector 
field From da;''/dA = ^ generates on M. a congruence of curves x^{X): thus A is the parameter 
along the congruence. Given a point p, this will always lie on one of these curves, and we can always 
take p to correspond to A = on this. The coordinates of a second point (7 at a parameter distance A 
from p on the same curve, will be given by 

i^'{X)^x^' + + ■■■ , (Al) 

where the x^ are the coordinates of p and the are those of g, approximated here at first order in A. Eq. 



(Al) is usually called an 'infinitesimal point transformation', or an 'active coordinate transformation' 
(see, e.g. Ref. ||, page 70; Ref. [|6|, page 49; cf. also Ref. [0, page 291, and [|9), Appendix C). At 
the same time, we may think that a new coordinate system has been introduced on A^, defined in 
such a way that the y-coordinates of the point q coincide with the x-coordinates of the point p; using 



(Al) it then follows from this definition that 

y^^(q) := x^^ip) = x^{q) - \Cix{p)) + ■■■ 

^x'^iq)-Xe{xiq)) + --- . (A2) 

In practice, we have in this way defined at every point a 'passive coordinate transformation' (i.e., just 
an ordinary relabeling of point's names), that at first order reads: 

?;^(A) =x'' -X^" + ■■■ . (A3) 

Suppose now that a tensor field has been given on A4; e.g., to fix ideas, consider the vector field Z 
with components in the x-coordinate system. In the same way that we defined a new coordinate 



system y'^ once a relation between points was assigned through (Al) by the action of ^, so we can now 
define a new vector field Z , with components in the x-coordinates, such that these components at 
the coordinate point x^{p) are equal to the components Z'^ the old vector Z has in the y-coordinates 
at the coordinate point y{q): 



Z^{x{p)) := Z'^{y{q)) = 



dy^ 



dx" 



Z%x{q)) . (A4) 

x(q) 
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The last equality in this equation is just the ordinary (passive) transformation between the components 
of Z in the two coordinate systems: we need it in order to relate Z and Z in a sing le s ystem (the 
a;- frame here), thus eventually obtaining a covariant relation. Indeed, substitution of (A3) into (A_4) 
and a first order expansion in A about a;(p) in the RHS gives 



Z^(A) = Z^ + A£^Z^ + ^ 
£^Z>^ :-ZX-C^' ' 



(A5) 
(A6) 



where, given that the point p is arbitrary, the dependence of all terms by x(p) has been omitted. The 
vector field Z is called the pull-back of Z, because is defined by dragging Z back from q to p, an 
operation that gives at p a new vector with components Z^ ^ given by (A4). In the particular case of 
the transformation ( [A^ ) this is the Lie-dragging. Now, having at the same point two vectors, these can 
be directly compared: at first order, Z{X) and Z are related by (A5), ([A^). In fact, in the limit A ^ 0, 
is this comparison that allows us to define the Lie derivative, with components (|A6|); Eq. (A15) below 
generalizes this to a generic tensor T. 

recalling it in some detail allows us to 



Although the story so far is a textbook one (cf. [ p9y3 
easily extend it to higher order. First, one has to realize that ( Al) is just the first order Taylor expansion 
about x{j)) of the solution of the ordinary differential equation dx^/dA = defining the congruence 
a;'^(A) associated with ^. The exact solution of this equation is the Taylor series (cf., e.g., [p8|, page 43) 



x^{q) = x^[j>) + K'^i.Av)) + y e^.r(x(p)) + ■ 



(A7) 



on using dx^ / d\ = d'^x'^/dA'^ = ^'^'C'^, etc. In practice, since p and q are arbitrary, we may simply 
write 



^^(A) = x^+Ae^ + ye^.r 

= exp[A£5]a;'' . 



(A8) 
(A9) 



The latter exponential notation is useful, in that it allows us to see the coordinate functions as the 
pull-backs of the functions x^ given by the exponential pull-back operator exp[A£^]. Furthermore, it is 



clearly seen by exp[(Ai + A2)=f4] — exp[Ai £^ exp[A2£^] that the point transformations (AS) form a one- 
parameter group of transformations. Using again the definition y^ijq) := x'^{p) for the y-coordinates, 
we get from (A7) 



j/'^(A)=a;^-Ae'' + ye',.r 



(AlO) 



on expanding all terms about x{q), eventually omitting again the x(q) dependence, since q is arbitrary. 
Finally, using ( A10| ) into (A4) and expanding all the terms about x{p), we get the a;-components 2'^ (A) 
of the pull-back Z{\), which reads [^ : 



Z^(A) = [exp[Ai'^] Z]^ 

= Zi' + \£.Z^ + — £f Z'^ • ■ • . 
Eq. ( |A4| ) is readily generalized to more general tensors than Z: we simply have to add to the RHS 



(All) 
(A12) 



of (A4) the right number of transformation matrices. Thus, the pull-back T of a tensor field T of type 
(p,q) is defined by having x-components given by 



f^--'^%,....,(x(p)) :=r'^--'^%,....,(y(g)) 



(A13) 



x{q) 



Using (|A10| ) as above then gives, omitting indices for brevity. 
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f(A) = T + X£^T+ — £| T 



(A14) 



To summarize, each of the diffeomorphisms forming a one-parameter group, as mathematicians cah 
the transformations generated by a vector field ^ and represented in coordinates by (A9), gives rise to 
a new field, the pull-back T(A), from any given tensor field T and for any given value of A. Thus r(A) 
and T may be compared at every point, which allows one to define the Lie derivative along ^ as the 
limit A ^ of the difference f(A) - T: 



At higher order we have 



d 

dA 



A=0 



T(A)=lim{(f(A)-r 



(A15) 



4t 



dA*= 



f(A) 



(A16) 



A=0 



On the other hand, the relation at each point between any tensor field T and its pull-back T(A) is 
expressed at the required order of accuracy by the Taylor expansion (A14). 



2. One-parameter families of transformations 



In order to proceed, considering more general point transformations than ( |A8| ) and more general 
Taylor expansions that (AM), some general remarks are in order. First, it should be noticed that the 
definition y^{q) '■= x^{p) for the y-coordinate system is completely general, given a first coordinate 
system (the x-frame here) and any suitable association between pairs o f po ints (more precisely, any 
diffeomorphism) , of which the one-parameter group of transformations (A8) is a part icular example. 
Second, the same generality is present in the definition of the pull-back, Eq. (A13), which is also 
independent from the specific type of transformation chosen. 

As we said in Section IL exact gauge transformations do not form a one-parameter group, but a 
one-parameter family ||l^ , [lq| . However the consequences of this fact show up only with non linearity, 
which is why at first order gauge transformations are approximated by (Al), (A3) (cf. [ p9| , ^ -p7[ ) . 
Therefore, having in mind a second order treatment of perturbations, the question we now have to 
deal with is twofold: i) which is the general form of families of transformations that depend on one 
parameter (one-parameter families of diffeomorphisms), but do not form a group; ii) which is the form 
of the Taylor expansion of the pull-back T(A) of a tensor T generated by one such one-parameter family 
of transformations. 

In (cf. also 1 40 1) we have shown that the action of any given one-parameter family of trans- 

formations can be represented by the successive action of one-parameter groups, in a fashion that, to 
order A^, reminds us the motion of the knight on the chess-board: 



(A) 



A2 



(1) 



CM 

?(2; 



(A17) 



A vector field ^(j,) is associated to the fc-th one-parameter group of transformati ons, w ith parameter A^ 
(we denote Ai — A). Similarly to the knight, the action of the transformation ( |A17| ) first moves from 
point p (with coordinates x^^) by an amount A along the integral curve of ^(^ij [i.e., according to Eq. 
(A8)]; then, it moves along the integral curve of ^(2) by an amount A2 = A^/2. At each fc-th higher 
order, a new vector field ^(k) is involved, generating a motion by A^ ~ \^ /k\. Thus, the action of a 
one-parameter family of transformation s is a pproximated, at order fc, by a 'knight transformation' of 
order k (see Theorem 2), of whic h ( A17 ) is the second order example. 

Given the 'knight transformation' (A17), we can now use it to define the y-coordinates, which will 
be given by 



A^ 
2 



{x{p)) 



(A18) 



Expanding the various quantities on the RHS around and omitting the x{q) dependence, (A18) 
becomes finally 



23 



(A19) 



Using again the case of the vector field Z as our paradigmatic example, we can now derive the pull- 
back Z{X) generated by a one-parameter family of transformations. Substituting (A19) into (|A^), 
and expanding again every term about x{p), we obtain the cc-components Z'^(X) of Z(X), which (after 
properly collecting terms) at second order read 



Z'^(A) = Z'^ + Ai;^,„ Z^ + ^ [£l^ + £^^^^ 



(A20) 



For a generic tensor T, again omitting indices for brevity, use of ( A19 ) into ( [A13| ) obviously gives 



T 



(A21) 



APPENDIX B: SECOND-ORDER PERTURBATIONS OF USEFUL QUANTITIES IN THE 

SYNCHRONOUS GAUGE 

In this appendix we report the expansions up to second order of a number of tensors, which have 
been used in deriving the results of Section IV. All calculations are performed in the synchronous and 
comoving gauge, assuming an Einstein-de Sitter background. No subscripts s on synchronous- gauge 
quantities will be used in this appendix. 

The covariant conformal spatial metric tensor is expanded as follows, 

7.,=%+7g' + ^7j^ (Bl) 
The corresponding contravariant metric takes the form 



= S'^ - 7(i)*^' - ^T^")*^' + 7(i)'Sr^' , (B2) 



(12) 

where the indices of the perturbations ' are raised by 5^^ . 
The extrinsic curvature tensor i?*^ up to second order reads 



l/^„(i)*'j_i^(2)*'_^(i)»fc^(i) 



^^=2 +2^r-7^^^^'^7in- (B3) 



The square root of the metric determinant is 



with inverse 



2' ' 4'* 8''*^ 4 
^Froni these quantities we can easily get the density contrast 



i7«^^7{;^ - ^^'-^iS - ^7^ ^^^0 + , (B6) 



having assumed as initial conditions 7Q^j = and i5q^'' = (i.e. 5q ~ Sq'^)- 
The Christoffel symbols up to second order read 
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■pi 

1 41, 



(2)i (2)i _ (2),z 

^ j ,A; fe ,j ^jk 



(l)e , (1) 



7 ^ I T J ,fe + 7 fe J 



7j.^'^ , (B7) 



from which, after a lengthy but straightforward calculation, the conformal Ricci tensor of the spatial 
hypersurface. 



1 



(l)ifc (1)*: , 

7 ' ^ j 



1 



+ 7(7^'^".+7?:L-vVf-7^^^: 



1 

+ 2 
+7 



' (IM /^y2 (1) , (1)^ _ (IK V ,^(1).' 

7 I ^ 7fej + 7 I jfe 7 fe I J ,£k j +1 I 7 J ,fe^ + 7 J 7 fc 7 j£ ,fe 



{i)ei (i),fe (i)« (i)fe 



(l)^fe (l),i 



7 j ,fc + 7 fc J 7jfc '1+7 ,fc7 ij' 7 ,fe7 j ,^ + ^7 ,j7 ^fc 



H — 7 
2 ' 



and its trace 



(l)m / . (1) 



7 ,,• + 7 



7 



'^ — 7 7fc ,^ +2!' ,^fe Ik ,e j^i 7jfc +7^,jfe '^ij,ek 

, mek ( , 3 (i)^^ (i),fe 1 (i)e,j (i)k 1 (i)i,^ (i)fc 



(B8) 



(B9) 



follow. 



APPENDIX C: SECOND-ORDER PERTURBATIONS GENERATED BY LINEAR TENSOR 

MODES 

In this appendix wo report the equations which allow to determine the second-order perturbations 
which arise due to the presence of tensor modes at the linear level: these are originated both by the 
coupling of scalar and tensor modes and by tensor-tensor mode couplings. All calculations are performed 
in the synchronous and comoving gauge, assuming an Einstein-de Sitter background. 

The equations which follow refer only to those parts of the second-order metric perturbations which 
involve tensor modes in the source terms (hence the subscript (t)). 

Raychaudhuri equation 

,(2) "2 (2) ' _ ^ (2) i,y2 (1)T _ 1 (l)Ti/ (l)T' Ay(l)Ty' (1)T 

-^X^^^^-V^xir + - ^'^''^) - ; (CI) 

energy constraint 



^fsit) 3^ <Ps(t) + ^2^s(t) 12^^W'' 



_ J_y(l)Tu' (l)T' _ Ay(l)Tu' , 1 (l)Tyv2y(l)T 

+ }L^imiJ,kAm _ J_ (l)T*j,fc (1)T _ J_ (l)Tzj (1)T _ (l)T*j (1)T , . f^^. 
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momentum constraint 



2^(2) ',1 (2)» '__ 



(l)Tifc ' (l)Tfe,/\ , 1 



evolution equation 



^<^^™(xS7'-xL:!7')4x^^^™,x!r'; (C3) 



yX<rvV + f X(^^™^,«^^ - yX«™^.. - yxS"^'" + yV^X^JV 

•T"^ / (l)Ti.ke , (l)Tfc£.i (l)Tfci/ (l)Tfc.i<?\ 

,(l)T«fe'^WT' 



(l)TjA; (l)Tfe,A /5 , \ 



, / (l)Ti.fc (l)TjA; 



-^x«™'xir'^^ + x(^ 



V X <y5,fc£'3 J I - Y^X 'f,khnd j 

Xl)Tkl( {1)T^ , (1)T,» _^(1)T» \ 

IX A; j£ 1^ X fej X j ,fe£ X fcf ,j I 

''":ex''rj', ■ (C4) 

reen method. The 



i^gX km " ] ^x 



1 (l)Tfcm I 



~'~gX /A km j ^X ^^A m ^k j ' 

f2) 

The Raychaudhuri equation can be easily solved for ^g^t-j by means of the Green method. The 

resulting expression has been given in Eq. ( 4.32| ) of the main text. By replacing it in the remaining 

(2) 

equations one can in principle obtain the traceless tensor Xs(t)ij by integration. 

APPENDIX D: TENSOR CONTRIBUTION TO THE SECOND-ORDER GAUGE 

TRANSFORMATION 



In this appendix we show how to compute the contribution from linear tensor modes to the per- 
turbed metric in the Poisson gauge by performing a gauge transformation from the synchronous gauge 
perturbed metric obtained from Appendix C. The equations for the gauge transformation parameters 
involved are obtained straightforwardly from Eqs. ( p3.12 ) to (5.14) 



(1)T/ , 4 (1)1 1 



12 ^^3'^ ^ 



3 



6 



y2^(2) _ ^2fl(2)/ 



2t 



(Dl) 



^ "(t)j 



3^ P(t),j Xs{t)z,j 



2^ ^ 



T- (1)T„2 ,7 T' (1)T ,jk , 2t (i)t 



3 



Xjk 'P/ 
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The energy constraint (C2) can be used to replace xfil'^^j term of '/'g'^t) and products of first-order 



quantities. 

On the other hand, we have from Eqs. (5.8) to (5.11) that the contributions to the perturbed metric 
that we are interested in are 



^p(t) 


(2)1 , 2 (2) 
= «(t) + 7"(t) : 




(D2) 




2t (i)t (2) , 


P{t).i ^ "(t)i ' 


(D3) 


'Pp(t) 


- '''s(t) 9 Xy 




(D4) 


(2) 
Ap(t) ij 


(2) , 2t / (i)T 







-^^.■xll^V^"^ + 2 + D.,/3;g) . (D5) 

We can obtain expressions in terms of the synchronous gauge perturbed metric as follows: 
Lapse perturbation 

Replacing the expression for a^^j* into Eq. ( pl)2| ) we obtain an expression in terms of ^g^j) , its derivatives 



and products of first-order quantities, that with the help of the Raychaudhuri equation (CI) can be 
written as 



2t' 



+ il,ky2 (1)T _ 3 (1)T. (l)T/ kj,t . 3 (i)T (i)T, kJ,^ _ ^ (1)T fej.*y2 (1)T 

+ ^ (2x«^ - V^X?" + 2xg^ XW" - 2V^x^^l7x(^)^ ^ - 2x^i^X«" • (D6) 
Shift perturbation 



^From Eq. ( p3| ) and using the momcntmn constraint (C3) and the Raychaudhuri equation (CI) we 
obtain 



(-4X?"V- - 2x£7'x'^^^ + ^XSJ'X^^ + ^'X^" 
(4xif V'*^^- - xif V^x*^)^ - 2x(^)^ V^xif 



-3xS7'x^^^" '^-^ + 2xg7'x^^)^ 



(D7) 



Spatial metric, trace 



^From Eq. (D4) and using the Raychaudhuri equation (CI) we obtain 
1 



-P(t) 



s(t) ' 8 



V2(;^(l)^ 



-I 
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Spatial metric, traceless part 



Replacing from Eqs. (Dl) the expressions for a^^? and d^^y. in Eq. ( [[)5| ) we obtain 



72y72,.(2) 



1 



(t) 



1 



p{t)ij 



s(t)y 



(t) k{i,j) ' 2^(*) 2 



4r 

T 



4r^ 



2t2 



4r^ 



3 

2t2 



(1)T ,ifc 1 (1)T ,lk 



3) 



(D8) 



(D9) 
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